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Представлена методика расчета квазистатического пространственного напряженно-деформированного состояния ведущего уст- 
ройства при метании стержня. Методика позволяет моделировать его поведение при интенсивном нагружении и в некоторой 
мере заменяет дорогостоящие экспериментальные исследования по отработке конструкции сборки, подбору материалов, вели- 
чин внешних нагрузок относительно дешевым и оперативным численным экспериментом. Задача упругопластического дефор- 
мирования решается с помощью вариационно-разностного метода. Физические соотношения принимаются согласно теории 
малых упругопластических деформаций. Г еометрические соотношения берутся в виде уравнений Коши. Физически нелинейная 
задача решается методом переменных параметров упругости. 


1. Введение 

Процесс высокоскоростного метания характе- 
ризуется высокой интенсивностью действующих 
нагрузок на метаемые элементы и вспомогатель- 
ные, т.н. ведущие устройства (ВУ), с помощью ко- 
торых осуществляется разгон сборки. Повышение 
скорости метания элементов приводит к увеличе- 
нию действующих на сборку нагрузок, нежелатель- 
ному уровню деформаций и напряжений и, как 
следствие, разрушению, как самой сборки, так и 
метательной установки (МУ). 

Вследствие того, что экспериментальные исследо- 
вания компоновочных схем подобных конструкций 
являются необходимым и обязательным этапом при 
их проектировании и отработке, и учитывая большие 
временные, финансовые и экономические затраты 
при проведении этих работ, естественным является 
стремление проектировщиков и конструкторов ис- 
пользовать математические методы механики дефор- 
мируемого твердого тела для моделирования напря- 
женно-деформированного состояния (НДС) элемен- 
тов конструкции при метании с целью сокращения 
экспериментальных испытаний. Поэтому наличие 
методик анализа НДС подобных конструкций, осно- 
ванных на современных численных методах, является 
обязательным при их проектировании. 

В работе анализ деформирования ВУ при мета- 
нии проводится в квазистатической постановке, 
обоснование которой представлено в [1]. Для рабо- 
ты с появляющимися согласно принципу Даламбе- 
ра [2] фиктивными массовыми силами вводится 
понятие коэффициента перегрузки (или просто пе- 
регрузки У), который определяется отношением 
ускорения тела а в данный момент к ускорению си- 
лы тяжести & N = а/щ. Получаемая при этом систе- 
ма нагрузок является самоуравновешенной, а из 
рассмотрения исключается постоянная составляю- 
щая осевых перемещений конструкции. 

Определение трехмерного упругопластического 
НДС проводится вариационно-разностным мето- 
дом (ВРМ) [3-5] с использованием выражений для 
представления пространственных производных че- 
рез интеграл по замкнутой поверхности. Реализация 


системы линейных алгебраических уравнений боль- 
шого порядка относительно искомых перемещений, 
к которой приводит использование необходимого 
условия экстремума сеточного аналога функциона- 
ла полной потенциальной энергии системы, прово- 
дится итерационным методом верхней релаксации с 
использованием оптимального коэффициента ре- 
лаксации. Физически нелинейная задача решается с 
помощью соотношений деформационной теории 
пластичности и метода переменных параметров уп- 
ругости. Нелинейная зависимость интенсивности 
напряжений от интенсивности деформаций ап- 
проксимируется ломаной двухзвенной линией. 

2. Расчет НДС секторов ВУ 

Рассматривается задача определения НДС состо- 
ящей из стержня и ВУ сборки (рис. 1) при метании. 
Причиной пространственной постановки задачи яв- 
ляется конструктивная трехмерность ВУ. На левом 
торце сборки приложено давление газа Р, посред- 
ством которого она разгоняется в трубе МУ. Задача 
решается в декартовой системе координат х, у, і. 



ВУ типа "катушка" состоит из трех одинаковых 
секторов РРЫНОМ (рис. 2, а), между гранями кото- 
рых находится герметик, предотвращающий прорыв 
газов. Сектора ВУ со стержнем крепятся по части 
контактной поверхности "стержень-сектор ВУ" ра- 
диусом К 2 , ограниченной осевыми координатами І н 
и І в , с помощью резьбового соединения, называе- 
мого гребёнкой, и которое здесь не рассматривается. 
На сделанном разрезе одного из секторов ВУ тем- 
ным цветом и величинами 2 Н и 2 В осевой координа- 
ты г обозначена зона гребёнки, включающая в себя 
часть поверхностей стержня и секторов. 
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Стержень, выполненный из твердого сплава, 
считается недеформируемым и его НДС не анали- 
зируется. Но при постановке задачи его масса, ко- 
торая равна 2 кг, учитывается при определении ко- 
эффициента перегрузки N осевых массовых сил и 
приложенного на гребёнке сдвигового напряжения 
СБВ, появляющегося при движении сборки в трубе 
МУ. Исследуемые сектора ВУ изготовлены из изот- 
ропного материала с меньшей плотностью, нежели 
плотность материала стержня. 




Рис. 2. Г еометрия секторов ВУ и вариант нанесения конеч- 
но-разностной сетки 


Геометрия секторов ВУ задается параметрами: 
В І - внешний радиус секторов; В 2 - внутренний ра- 
диус секторов (радиус стержня); К 3 - радиус цили- 


ндрической части секторов; Я - длина секторов ВУ; 
Ь ь Ь 2 , Я 3 , Х 4 , Ц - параметры, характеризующие ге- 
ометрическую форму секторов. Описанная модель- 
ная задача реализуется для следующих значений ге- 
ометрических параметров: = 5,0 см; К 2 = 1,0 см; 
К, = 3,5 см; Ь = 15,0 см; Я, = Я 2 = Я 4 = Я 5 = 2,0 см; 
Я, = 7,0 см; 2 Х = 2,0 см, І 2 = 4,0 см, 7 3 = 11,0 см, 
7 4 = 13,0 см, І н = 5,5 см, 7ц =7,5 см 2 В = 9,5 см. Дли- 
на гребёнки вдоль оси 0<- составляет 4,0 см. 

Ввиду равноправности трех секторов при нагру- 
жении, а также симметрии параметров НДС секто- 
ра относительно плоскости Охі, при численном ана- 
лизе (как и представлении полученных результатов) 
рассматривается лишь половина МИНС сектора, 
величина которого задается углом <р~6(У (рис. 2, а). 

Задача решается при следующих статических и 
геометрических граничных условиях для половины 
сектора: 

нижний торец сектора г = 0: 

ст«=<7* = 0, а, = Р; (1) 

верхний торец і = Е\ 

о хі = <у ух. = оу = 0; (2) 

контактная поверхность "стержень-сектор ВУ", 
кроме гребёнки, (жесткая стенка со скольжением): 
и = ѵ = О, 

О" л’ Т1\ 1 СТ уу П у О, ((Г гх 0) 

гребёнка (жесткая стенка без скольжения): 
и = ѵ = 0, 

а хх п х + ОухПу = СОВ- (огх = СОВ) 
грань МИ (у = 0): 

ѵ = 0 , а ух = о ух = 0 ; 

контактные поверхности "труба-сектор ВУ" АВ и 
ЯМ (рис. 2, б) (жесткая стенка со скольжением): 
и = ѵ = 0, 

О" л’ П х + ОухПу 0, (ст хх 0) 

свободные поверхности ВС и БЕ (равенство нулю 
проекций на оси координат действующего на нак- 
лонной площадке полного напряжения): 

СУхПх "Ь СТ ху ТІ у Т" Схх/Іх 0, 

ОухПх Т" СТ уПу “Ь СТ ух Ну — О, 

(7хх п х + <7-уПу + (7х п х = О, 

свободная поверхность СО, задаваемая радиусом Я 3 - 

а х п1 + ОуП у + 2охуП х Пу = 0, (<т, = 0) 

(Оу-СТх)п х Пу +СТ Х у(Пх -Пу) = 0, (,0 г<р = 0) 

СТ хх Пх СТ ух Пу 0, (сТгх 0) 

свободная поверхность грани НС, для которой 
ср=6 0° (рис. 2 , а) - 

а х п 2 у + ступ 2 - 2а хуПхПу = 0, (аѵ = 0) (3) 

(ст у СТ х ) ПхПу + СТ Х у (и х Сіу ) 0, (ст ГІ р 0 ) 

СТ ух И х СТ у, Пу — О, (.СТ х(р 0) 
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где а„ а у , а ѵ ст ху , а уѵ а 7Х , о ух , от,., а п - компоненты 
тензора напряжений; и, ц \ѵ - компоненты вектора 
перемещений в декартовой системе координат; п х , 
п у , п ѵ - направляющие косинусы; сг г , а ѵ , о ѵ ст„ ; „ а п , 
су.„ ст (? ,., ст-р„ ст,.. - компоненты тензора напряжений в 
цилиндрической системе координат г, ср, і. 

Граничные условия на боковых поверхностях 
сектора ВУ в декартовых координатах получены с 
помощью формул преобразования компонентов 
тензора напряжения для цилиндрической системы 
координат. Поэтому для наглядности рядом с неко- 
торыми статическими граничными условиями в 
скобках содержатся записи этих же условий в ци- 
линдрической системе координат. На осевые пере- 
мещения \ѵ не накладывается ограничений. 

Представленные ниже результаты получены на 
конечно -разностной сетке (Ы]хк) = (15x17x66) при 
реализации системы примерно 50500 линейных 
алгебраических уравнений методом верхней релак- 
сации с выбором оптимального коэффициента ре- 
лаксации щ ор1 . Здесь і - количество узлов сетки по 
координате х(г) (по толщине секторов),) - количе- 
ство узлов сетки по координате у(ср ) (по окружнос- 
ти секторов), к - количество узлов сетки по коорди- 
нате г (по высоте секторов). В качестве критерия 
окончания итерационного процесса брались усло- 
вия выполнения теоремы Клапейрона или равен- 
ства искомых величин перемещений в двух сосед- 
них итерациях и упругих соседних задачах с заранее 
заданной точностью. 

Материал секторов - дюралюминий с плот- 
ностью р = 0,0027 кг/см 3 и значениями упругих 
констант: коэффициент Пуассона р =0,3; модуль 
упругости Е =70000 МПа; модуль упрочнения 
Е і = Е/5; деформация начала текучести е 5 = 0,003. 
Масса ВУ для заданной геометрии и выбранного 
материала равна 2,066 кг. 

Величины действующих на сектора перегрузки 
и сдвиговых напряжений на гребёнке определяются 
из условия равенства нулю проекций всех внешних 
сил на ось 0г и для заданной величины давления 
Р = -250 МПа составили N ~ -48254,8 и 
СБВ ~ -352,8 МПа соответственно. Направление 
действия нагрузок N и СРВ противоположно нап- 
равлению движения сборки (направлению оси 0ф. 

Распределение осевых напряжений а ѵ а также 
интенсивности напряжений ст, и сдвиговых напря- 
жений а п по длине секторов (от значения осевой 
координаты г = 0 до і = Е) вдоль плоскости сим- 
метрии Охг {(р = 0, грань МУ на рис. 2, а) у их внут- 
ренней поверхности г = В 2 представлены в таблице 
в строках п = 1, 2, 3. 

В силу характера действующих внешних нагру- 
зок осевые напряжения изначально являются оп- 
ределяющими в НДС метаемой сборки. Практи- 
чески во всем объеме сектора ВУ эти напряжения 
сжимающие. Исключением является область на 
контактной поверхности "стержень-сектор ВУ" в 
районе гребёнки, где метаемый стержень крепится 
к секторам ВУ. Наибольшее значение напряжения 


сг, = -795 МПа достигается непосредственно под 
гребёнкой. В нижней части гребёнки осевые нап- 
ряжения су. = -758 МПа, будучи сжимающими, ме- 
няют знак и становятся растягивающими, достигая 
величины 385 МПа в самом ее верху. Помещенные 
в упомянутых строках таблицы другие напряжения 
также имеют максимальные значения на гребёнке. 


Таблица. Распределение напряжений в ячейках у внутрен- 
ней поверхности сектора на плоскости симметрии 
Охг ВДОЛЬ ОСИ 02 


п 

Напря- 

жения, 

МПа 

Осевые координаты і, см 



0,050 

4,875 

5,375 

5,625 

6,125 

7,375 

9,375 

10,10 

14,80 

т 

л- 

-250 

-591 

-795 

-758 

-431 

-118 

385 

94,6 

-2,8 

2 

<7, 

251 

472 

665 

733 

646 

558 

505 

231 

27,5 

3 

а* 

-1,0 

-14,9 

-85,6 

-216 

-293 

-310 

-228 

-27,6 

1,0 

4 


-250 

-398 

-391 

-384 

-363 

-294 

-204 

-187 

-5,5 

5 

О, 

245 

282 

297 

302 

304 

277 

177 

148 

34,6 


Полужирным шрифтом выделена зона гребенки и относящиеся 
к ней результаты 


Окружные напряжения сг р , также как и радиаль- 
ные напряжения ст п являются сжимающими прак- 
тически во всем объеме секторов. От нижнего торца 
постепенно, по мере увеличения осевой координаты 
I, напряжения ст,, возрастают по абсолютной вели- 
чине, достигнув на внутренней поверхности г = К 2 
секторов в ячейке с величиной угла ср ~ 0 и осевой 
координаты г = 5,375 см, т.е. прямо под гребёнкой, 
наибольшего значения су 0 = -317 МПа. Далее (по 
высоте сектора) окружные напряжения достаточно 
быстро уменьшаются. 

О характере поведения осевых су. {а) и окруж- 
ных су,, (б) напряжений по толщине секторов ВУ 
можно судить по представленным на рис. 3 резуль- 
татам для нижней части гребёнки г = 5,625 см. Нап- 
ряжения довольно резко уменьшаются по величине 
от внутренней поверхности секторов к их внешней 
поверхности г = В,. Окружные напряжения, кроме 
того, уменьшаются от грани МЫ, для которой угол 
ср = 0, до нулевых значений на свободной поверх- 
ности - грани НО, как того требуют заданные гра- 
ничные условия (3). По поводу выполнения стати- 
ческих граничных следует отметить, что одним из 
достоинств ВРМ является автоматическое выпол- 
нение этих условий (также как и выполнение урав- 
нений равновесия) при минимизации функциона- 
ла полной потенциальной энергии системы. 

Результаты на рис. 3 были построены следующим 
образом. В плоскости Оху для фиксированного значе- 
ния координаты г на перпендикулярах из середин 
ячеек откладывались значения напряжений, которые 
для наглядности соединялись в двух направлениях 
линиями и образовывали некоторую поверхность. 

В зону пластических деформаций попадает око- 
ло 51 % объема секторов, включая примерно две 
третьих части (по высоте) гребёнки. На всю толщи- 
ну сектора материал деформируется пластически 
от координаты г ~ 2,7 см до г ~ 8,35 см. Значения 
интенсивности деформаций е, увеличиваются при 
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Рис. 3 . Распределение а) осевых <т 2 и б) окружных сг ѵ напряжений в координатной плоскости Оху сектора ВУ в нижней части 
гребёнки і = 5,625 см 


приближении к гребёнке, как от торца г = О, так и 
по радиусу от внешней поверхности сектора. 

На рис. 4 кривой 1 представлено распределение 
осевых напряжений ст, на внутренней поверхности 
секторов вдоль оси г для значения угла ср = 0. На гре- 
бёнке (она обозначена штриховкой) имеет место 
большая неравномерность распределения внутрен- 
них напряжений. При движении сборки в трубе МУ в 
нижней части гребёнки возникают большие сжимаю- 
щие осевые напряжения, приложенные к нижним 
виткам резьбы. В верхней же ее части - несколько 
меньшие по величине растягивающие напряжения, 
которые приложены к верхним виткам. Амплитуда 
изменения этих напряжений вдоль гребёнки длиною 
4,0 см составляет Асг, = 1180 МПа. Наибольшие зна- 
чения величин интенсивности деформаций секторов 
находятся на этом участке их внутренней поверхнос- 
ти. Максимальное значение параметра (е,) шах ~ 0,0404 
имеет место в нижней части гребёнки для значения 
осевой координаты г = 5,625 см и угла <р ~ 54,4° вбли- 
зи грани НО. Поэтому, учитывая резкое изменение 
осевого напряжения на гребёнке, и то, что остальные 
напряжения достигают наибольших значений также 
на гребёнке, можно сделать вывод о том, что этот 
участок внутренней поверхности секторов является 
наиболее опасным с точки зрения нарушения цело- 
стности ВУ и всей сборки при метании. 

Описанная неравномерность распределения нап- 
ряжений по величине и направлению характерна для 
небольшого объема сектора, находящегося у гребёнки, 
и носит локальный характер. Из представленных на 
рис. 3 результатов видно, что по мере удаления от гре- 
бёнки в направлении внешней поверхности г = и 
грани НС напряжения достаточно быстро уменьшают- 
ся по величине. Так, при удалении от внутренней по- 


верхности по радиусу на 0,44 см (толщина сектора этой 
его цилиндрической части равна 2,5 см), наибольшие 
значения напряжения су 1 уменьшаются до величин 
- 565 МПа и 97 МПа в нижней и верхней частях гре- 
бёнки соответственно, а амплитуда - до величины 
Асг, = 660 МПа. Выравнивание значений парамет- 
ров НДС секторов наблюдается также и при удале- 
нии от гребёнки в сторону нижнего и, особенно, 
верхнего торца. 

В строках п = 4, 5 таблицы и кривой 2 на рис. 4 
представлены результаты, полученные при отсутствии 
стержня. В этом случае нагрузка на гребёнке СОВ = 0, 
а перегрузка /Ѵ~ -91167,1. Напряжение <т. без переме- 
ны знака согласно граничным условиям (1) и (2) плав- 
но меняет свою величину по длине секторов от нижне- 
го торца г =0, где ст, = -250 МПа, до практически нуле- 
вого значения на верхнем торце і = Ь. Аналогично по- 
ведение интенсивности напряжений о;. Сдвиговые 
напряжения <з к на поверхности практически отсут- 
ствуют. В зону пластических деформаций попадает 
41 % объема секторов. Величина (е,) тах ~ 0,0142 при- 
мерно в 3 раза меньше, нежели при наличии стержня, 
и достигается в месте сопряжения конической поверх- 
ности ВС и цилиндрической поверхности СП у внеш- 
ней поверхности сектора для ср ~ 30°. 

Таким образом, полученные результаты позволя- 
ют сформулировать следующий вывод. Причиной об- 
разования сложной упругопластической картины де- 
формирования материала с изменяющими на неболь- 
шом пространстве у гребёнки свои знаки напряжени- 
ями и деформациями и приводящими к образованию 
зон пластических деформаций сложной конфигура- 
ции, распространяющихся на всю толщину секторов 
ВУ, является приложенное на гребёнке внешнее сдви- 
говое напряжение СБВ от метаемого стержня. 
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Для изготовления секторов ВУ могут приме- 
няться конструкционные материалы с меньшими 
значениями модуля упрочнения Е { . На рис. 4 кри- 
вая 3 получена для Е ] = Е/25. Величина амплитуды 
на гребёнке Дет. = 1380 МПа. По сравнению с рас- 
смотренным выше вариантом задачи величины 
напряжений изменились незначительно, но дефор- 
мации выросли примерно в 4. ..4, 7 раза. Пластичес- 
ки деформируется 61,6 % объема секторов ВУ. На- 
ибольшая интенсивность деформаций увеличилась 
почти в 4,5 раза - (е,) тах « 0,1775, что совершенно 
недопустимо для метаемой сборки. 
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Рис. 4. Распределение осевых напряжений оу на внутренней 
поверхности секторов вдоль продольной оси і для 
ср = 0. Штриховкой обозначена гребёнка. 1) основной 
вариант задачи (дюралюминиевые сектора); 2) ос- 
новной вариант при отсутствии метаемого стержня; 
3) основной вариант при Т, = Е/25; 4) сплав В95 [7] 

Одним из путей понижения уровня напряжений 
и деформаций в секторах ВУ при метании является 
использование материалов с более высокими зна- 
чениями прочностных характеристик, которые 
"... должны: 1) иметь невысокие плотности и 2) вы- 
держивать большие деформации не разрушаясь" [6]. 

Был проведен расчет НДС секторов, выполнен- 
ных из сплава В95 [7], для которого Е х = Е/ 50, 
е. = 0,0085. Вследствие того, что величины плотности 
материала и давления газа остались прежними, зна- 
чения перегрузки и напряжения на гребёнке также не 
изменились. Распределение осевого напряжения а, 
по длине гребёнки представлено на рис. 4 кривой 4. 
Величина осевых напряжений на гребёнке уменьши- 
лась, их амплитуда Да, = 920 МПа. Максимальное 
значение интенсивности деформаций достигается в 
том же месте, что и для дюралюминиевых секторов, а 
его величина (е,) тах ~ 0,0142 стала меньше почти в 
3 раза. Зона пластических деформаций толщиною 


около 0,2 см располагается на внутренней поверхнос- 
ти секторов от значения координаты г = 5,375 см до 
г = 6,375 см. Всего 0,57 % объема секторов деформи- 
руются пластически. Таким образом, использование 
сплава В95 позволяет уменьшить на гребёнке значе- 
ния напряжений и особенно деформаций. 

Одновременно с применением новых материа- 
лов для конструирования секторов следует прово- 
дить исследование влияния их геометрии на НДС 
ВУ. Так, для понижения уровня напряжений и де- 
формаций в секторах ВУ при метании стержня мож- 
но предложить достаточно очевидное техническое 
решение - увеличить площадь поверхности контак- 
та метаемого стержня и секторов и (или) толщину их 
цилиндрической части. В последнем случае следует 
ожидать уменьшения величины осевой перегрузки, 
т.е. уменьшения скорости метания сборки. 

3. Выводы 

Созданы физическая и математическая модели 
для расчета трехмерного квазистатического упру- 
гопластического напряженно-деформированного 
состояния секторов ведущего устройства при мета- 
нии стержня при нагружении их комбинацией пове- 
рхностных и массовых сил большой интенсивности. 

Реализуемая в разработанном варианте вариаци- 
онно-разностного метода методика решения прост- 
ранственных упругопластических задач расчета ве- 
дущих устройств с использованием теории малых 
упругопластических деформаций и метода перемен- 
ных параметров упругости дает возможность оцени- 
вать напряженно-деформированное состояние для 
случая диаграммы растяжения материала, достаточ- 
но близкой к пластическому течению. 

Численный анализ деформирования ведущего 
устройства при метании выявил зону гребёнки, как 
наиболее нагруженный локальный участок секторов, 
где отмечается значительная неравномерность расп- 
ределения напряжений и деформаций, способная 
привести к разрушению сборки. Причиной такого 
распределения параметров напряженно-деформи- 
рованного состояния на гребёнки является нагрузка 
от метаемого стержня. Использование для секторов 
сплава В95 с более высокими прочностными харак- 
теристиками привело к значительному снижению 
уровня напряжений и деформаций на гребёнке. 

Следует отметить, что созданная методика поз- 
воляет получать информацию обо всех параметрах 
пространственного напряженно-деформированно- 
го состояния с точностью до ячейки конечно-раз- 
ностной сетки при проектировании и оценке проч- 
ности не только рассмотренного варианта сборки, 
но и других конструкций для широкого спектра 
внешних нагрузок. 
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